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Abstract 

Preparing many body entangled states efficiently using available interactions is a challenging 
task. One solution may be to couple a system collectively with a probe that leaves residual entan- 
glement in the system. We investigate the entanglement produced between two possibly distant 
qubits 1 and 2 that interact locally with a third qubit 3 under unitary evolution generated by pair- 
wise Hamiltonians. For the case where the Hamiltonians commute, relevant to certain quantum 
nondemolition measurements, the entanglement between qubits 1 and 2 is calculated explicitly for 
several classes of initial states and compared with the case of noncommuting interaction Hamilto- 
nians. This analysis can be helpful to identify preferable physical system interactions for entangled 
state synthesis. 
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I. INTRODUCTION 



Recent advances in the preparation and manipulation of many body entangled states 
has promoted keen interest in characterizing entanglement as a resource for quantum in- 
formation processing with applications to such areas as precision measurement, quantum 
communications, and computation Given a many body system Hamiltonian it is possi- 
ble to construct designer entangled states using only single body and at least one two body 
interaction, however, it is difficult in general to find the optimal sequence of control opera- 
tions to prepare designer entangled states. Often the number of control operations needed 
scales poorly with the number of bodies, or interaction with the environment introduces 
significant decoherence during the preparation. One candidate for efficient production of 
multiparticle entanglement is to collectively interact members of a system with a probe that 
can be measured after the interaction. This has been proposed in the context of quantum 
nondemolition (QND) measurements wherein a probe p interacts with the many particle 
system s via an interaction Hamiltonian H ps . If an observable M s of the system commutes 
with Hp S and is a constant of motion in the absence of the interaction, then measurement 
of an observable M p that does not commute with with H ps constitutes a QND measurement 
of M s . Such a technique can be used to prepare many body entangled states of the system. 

As an example, the interaction of a single frequency mode (two polarization mode) laser 
field in a coherent state with a collection of two level atoms can be modeled in a truncated 
Fock space as a 2n-dimensional system interacting with a collection of qubits. A nonde- 
molition measurement of the z component of collective spin of the atoms, S z = J^i^l, is 
made by interacting the field with the atoms collectively via a pairwise commuting interac- 
tion H = gJ2iJz^l, where J z = l/2(a+a + — a!a_) is the angular momentum component 
along z in the Schwinger representation where a± are the annihilation operators for each 
polarization mode. A subsequent measurement of a canonically conjugate variable of the 
field such as J x = l/2(a+a_ + a)_a + ) completes the QND measurement. Such an interac- 
tion can be realized as the polarization dependent AC stark shift on the atoms ||. Under 
unitary evolution the atoms become entangled with the field and by measuring the field it 
has experimentally been shown that one can prepare entangled atomic spin squeezed states 
||. To see this consider 2 two level atoms, 1 and 2, that interact with a photonic qubit 3 
storing information in the polarization degrees of freedom. Let 3 be prepared in the state 
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I J x = J = 1) = l/y/2(a\ + al)\Q) where |0) is the vacuum state, and the two qubits be 
prepared in a product state with both spins oriented along x, or written in the coupled basis 

\s. = i) = \ U\W 

= i(|l,-l) + V2|l,0) + |l,l). 
Under evolution by the Hamiltonian H = gJ z ® S z the state evolves to, 

U(t)\S x = l)\J x = l) = I(cos(f)(|l,-l) + |l,l))|J, = i) + v / 2|l,0)|J, = i) 

+ishi(f)(|l I -l>-|l,l»|J, = -i». 
After interaction times satisfying gt = ir(2m + 1), for m any natural number, a measurement 
of 3 along \ J X = 1) yields the conditional maximally entangled state for 1 and 2: |1,0) = 

|V + ) = i/\/2(iT,>iU + |T z >i| z ». 

It would be useful to know what kind of entanglement is produced within the collective 
sample before a measurement is made on the probe. Such knowledge could indicate what 
type of generalized measurement is optimal to prepare the desired state of the system. 
It is generally quite difficult to characterize the dynamics of entanglement between many 
bodies. Indeed, even a tractable quantification of entanglement between one qubit and an 
n dimensional system for n > 2 has yet to be found in the general ||. We study the 
simplest model possible namely the dynamics of entanglement between two qubits 1 and 2 
that interact pairwise with a third qubit 3. We show that there are some easily identified 
constraints to the dynamics of entanglement starting from various initial pure states of the 
three body system in the case that the interaction Hamiltonians Hi 3 and if 23 commute. This 
is the regime encountered in QND measurements of the type outlined above. In the case of 
noncommuting Hamiltonians the dynamics are less restrictive but can be characterized by 
certain symmetries of the interaction. 



II. COMMUTING INTERACTIONS 



We begin with a system of two qubits 1 and 2 interacting with a third qubit 3 via pairwise 
interaction Hamiltonians: 

333 

#13 = E ® ^ + E a >l ® 13 + E "fc 11 ® 4, 

iJ 3 =1 ^ T (2-1) 

#23 = E / V* ® °j + E ® 13 + E ft? ® °l 

i,j=l fc=l fe=l 
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where ex™ denotes the Pauli operator ct( 12j 3) = 0( x ,y,z)i acting on the space of qubit m. If 
[#13,#23] = 0, then it is straightforward to show that the Hamiltonians must assume the 
following form: 

#13 = \u\<?l® erf + \a'W\ ® l 3 + a"l 1 ® af t s 

#23 = Wl ® o) + \0'\al ® l 3 + /?"1 2 (g) a 3 , 
where a? 1 = 1 /|tI (Ti* 7 ! 1 + + 73°"™) • The second two terms in each of if 13 and if 23 

commute with all other terms and have the effect of local operations acting after the evolution 
generated by the entangling Hamiltonians: H[ 3 = \a\a\ Cg) crj, H' 2Z = |/3|<x| ® a 3 . Local 
unitary operations leave the entanglement between qubits 1 and 2 invariant and as such 
we consider only the reduced state of qubits 1 and 2 after joint evolution by the unitary 
U(t) = e - i(H w+^2 3 )*. 

An arbitrary pure state of qubits 1,2 and 3 can be written in the eigenbasis of o"i®cr|(8>cr 3 
which is the simultaneous eigenbasis of ^13,^23, as 

l*) = E a ±±±l ± )s® l%®l±)i = E a ±±±l ±±± )- ( 2 - 3 ) 

After evolution for a time t, the state is compactly written 

U(t)\V) = E « mi , m2 , m3 e- l(mim3|a1+m2m3|/3l) Vi^2m 3 ), (2.4) 

mi,m2,m3e{— 1,+1} 

and the reduced state of qubits 1 and 2 is 
Pn(t) = Tr 3 [U(t)\*)(y\W(t)} 

all me{-l.+l} 

(2.5) 

Wootters has given an explicit formula to calculate the entanglement of formation £(p) 
for an arbitrary state of two qubits in terms of the tangle defined as 712 = [max{\i — A 2 — A3 — 
A 4 }] 2 where the Aj are the square roots of the eigenvalues, in decreasing order, of P12P12, where 
P12 is the spin-flipped version of pi 2 , p 12 = ay®a y p\ 2 ay®a y . The entanglement of formation is 
a monotonically increasing function of the tangle given by S(p) = ft,(l/2 + l/2\/I — r) where 
h(x) is the binary entropy function. A normalized pure state of three qubits will generally 
have 15 independent real parameters (5 of which are invariant under local unitaries ||) 
making it cumbersome to calculate this quantity. We focus on simply determining how the 
entanglement between 1 and 2 changes after joint evolution for different classes of initial 
states. 
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A. Fully separable initial states 



Consider the case when the initial state is pure and fully separable, 

= Ri ®R 2 ®R 3 \+++), (2.6) 

k 

where the local unitary operators are written R k = e llk °% . Quite generally, 



(2.7) 



Rs\+) = c\+) + d\-), 
RsH = d*\+)-c*\-), 

where |c| 2 + |<i| 2 = 1, and we consider the action of U(t) on the states R± <g> R 2 \++±) 
separately. Because [Hi 3 ,H 23 ] = 0, we write U(t) = e -lH i3*e -lfl 23* = Ui 3 (t)U 23 (t) . Using 
the fact that the operators in the generating Hamiltonians are unitary and hermitian the 
unitaries can be expanded as 

U 13 (t) = cos(|a|t)l -ism{\a\t)a l d ®a] 
U 23 (t) = cos(0\t)l - ism(0\t)al® a*, 

and similarly, R k = cos^jk) — ^ sm (7fc)°"^ fc - The evolution can be divided into two pieces, 

Ui 3 (t)R 1 \++±) = (cos(\a\t) cos(7i)l - icos(\a\t) sin(7i)cr^ 

-ism(\a\t) cos(7i)al (g) a] - sin(|a|t) sin^cricrl <g> cr|)|++±) 
= e ^l 5 l*4 J R 1 |++±) 
= V 1± (t)\++±), 

U 23 (t)R 2 \++±) = e T¥ HR 2 \ ++± ) 
= V 2± (t)\++±). 

(2.9) 

The evolved state of the three qubit system is generally entangled, 

U(t)\V) = cV 1+ (t) ®V 2+ (t)\+++) + dV 1 _(t) ®V 2 _(t)\++-). (2.10) 
Tracing over the third qubit, the reduced state of qubits 1 and 2 is 



Pn(t) = \c\ 2 v 1+ (t)®v 2+ (t)\++)(++\v 1+ (ty®v 2+ (ty 

+ \d\ 2 V 1 _(t) <g> V 2 _(t)\++)(++\V 1 _(ty ®V 2 _(rf; 
a convex sum over separable states and therefore seperable, £(pi 2 (t)) =£(pi 2 (0)) = 0. 



(2.11) 
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The mapping on the reduced state p\ 2 = Tr 3 (\^/)(^/\) can also be written 

Pn(t) = 5>i ± (£)4 ® V 2± (t)RU( ± |^)(^|±) 3 yi ± (t) t i?i ® ^(t)^, (2.12) 

where i?i and _R 2 transform the basis states given in the expansion of |\&) to the eigenbasis 
of the interaction. In the case that the state of qubits 1 and 2 is initially uncorrelated with 
qubit 3, i.e. \ty) = \x)i2 ® |0)3, then it is convenient to express the evolution of the reduced 
state p\ 2 {t) in Krauss form J!]] 

p 12 (t) = -£A ± (t)\ X }(x\A ± (t)\ (2.13) 

where the Krauss operators are defined 

A ± (t) = 3 < ± |*7(t)|*) 3 = 3< ± l^s^Jt)^ <g> V 2± (t)4, (2.14) 

satisfying the trace perserving condition: X) A±(i)M.±(£) = 1. One could also find a Krauss 
operator expansion for the case where the initial state of the joint system of qubits 1,2 
and qubit 3 are entangled. This would require finding the unitary that maps an initial 
unentangled state \x)i2 ® \</>)3 to the entangled one. 



B. Entangled initial states 

Here we consider the action of the unitary evolution on initially entangled states from 
various classes. 



1. Bipartite entanglement between 1 and 2 
Let the initial state be 

l*> = |X>12|0>3 

= (a|O) 1 |O) 2 + 6|l) 1 |l) 2 )|0) 3 , 

where we have expanding the entangled state \x) between qubits 1 and 2 in a Schmidt 
basis with real parameters a, b. This state has an initial tangle between qubits 1 and 2 
Tradx)) = 4(a6) 2 ). Using \<f>) = Rs\+) with R 3 given in Eq. |2~7| , the state |\&) can be written 
as a linear combination of eigenstates of 

= cR 1 ® R 2 (a\++) + b\ — »i2|+> 3 + dR 1 ® R 2 (a\++) + b\ — })i 2 |-} 3 - (2.16) 



The Krauss operators are therefore, 

A+{t) = cV 1+ {t)R\®V 2+ {t)R\ 
A_{t) = dV l _{t)R\®V 2 _{t)Rl 

with Vj ± given by Eq. |2]9] with the appropriate rotation operators Rk for the state The 
reduced state is 

Pl2 (t) = Y,A ± (t)\x)(x\A ± (t)l (2.18) 

By the convexity of the entanglement of formation |3J, £ (a) > J^iPi^i^i) where the state Oi 
is related to a by local operations and J2Pi = 1- Therefore, £{p\ 2 {t)) <£(pi2(0)). 



2. Bipartite entanglement between 2 and 3 

Consider the following initial state which is seperable over qubits 1 and 2 

l*> = \Mx)* (219) 

= |0>i(a|O) 2 |O) 3 + 6|O) 2 |O) 3 ). 

As before, we have expanded the entangled state between qubits 2 and 3 in a Schmidt 
basis with real parameters a, b. This state is equivalent (up to local rotations) to a linear 
combination of eigenstates of a\ (g> cr| (g> erf as 

|tf ) = aR x ®R 2 ® R 3 \+++) + 6i2i ® #2 ® #s| + — )• (2.20) 

Here a Krauss operator expansion is not particular helpful because the two qubit system 
1,2 is initially entangled with qubit 3. However, the evolved three body state is easily found 
using the expansion for R 3 , 

U(t)\V) = V l+ (t) ®V 2+ (t)(ac\++) + bd*\+-)) 12 \+) 3 

+V l _(t) <g> V 2 _(t)(ad\++) - bc*\+-)) 12 \-) 3 (2.21) 
= Ep±V 1± {t)^V 2± (t)\+) 1 \u ± ) 2 \±) 3 , 



where 



and 



p+ = ^J\ ac \ 2 + \bd*\ 2 ,p_ = ^\ad\ 2 + \bc*\ 2 , (2.22) 



\u+) = l/p+(ac\+) +bd*\-)), 
|u_) = l/p-(ad\+)-bc*\-)), 



(2.23) 
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are the normalized states for qubit 2. The reduced state is then, 

Pu(t) = £*4ti±(*) ® V 2± (t)\+u ± )(+u±\V 1± (ty ® V 2± (t)t ; (2.24) 

a convex combination of seperable states which is therefore seperable: £(pi2{t)) =£(pi2(0)) 
= 0. This argument is completely symmetric under interchange of 1 and 2 for the case of 
initial bipartite entanglement between qubits 1 and 3. 



3. Three party entangled states 

The third set of initial pure states to consider is states that cannot be seperated over 
any set less than order three. There are several classifications of three party entangled 
states, but we focus on a division imparted by a physical quantity introduced by Coffman 
et.al. || known as the residual tangle r 123 . This quantity can be thought of as the amount 
of entanglement between 1 and the joint system 2,3 that cannot be accounted for by the 
entanglements of 1 with 2 and 3 seperately. For qubits such a quantity is naturally symmetric 
with respect to interchange of particles and is shown to be equal to 

T 123 = 2(\\ 2 \l 2 + \?\l 3 ), (2.25) 

where and are the square roots of the two eigenvalues of pijpij- A notable result 
of Diir et.al. Jl0| shows that states with zero residual tangle cannot be converted to GHZ 
states (states with T123 = 1) under stochastic local operations and classical communication 
(SLOCC). They shown that all members of the set of zero residual tangle (ZRT) states can 
be written, 

|^« T ) = a |000) + 6|001) + c|010) + d|100). (2.26) 
Brun and Cohen have constructed a distillation protocol using local POVM's that 



evolves all ZRT states to a boundary of the set defined triple states 

\y tr ) = /|001) + 0|O1O) + h\100). (2.27) 



This is one class of initial states we study under evolution by Eq. |2~4] . The second class we 
consider is the set of generalized GHZ states, 



m GHZ ) = a |000) + 6|111), (2.28) 
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which have residual tangle r 12 3 = 4a6 and equal the standard GHZ state when a = b = 1/ y/2. 
While these are not the most general states with nonzero residual tangle, it is the class 
most easily distilled to a GHZ state [0]]. The states |\l/* r ) and \ty GHZ ) are inconvertible 
under SLOCC. 

a. Generalized GHZ states 
Let 

I^ghz^ = a |ooo) + 6|111) 

^2 . 29 J 

= a J Rii?2#3(a|+++) +6| ))• 

Tracing over any one of the qubits leaves a mixed state for the other two with no entangle- 
ment. Under the action of R3 (from Eq. 16), 

\m GHZ ) = R 1 R 2 a C \+++) + ad\++-) + bd*\ +) - bc*\ ). (2.30) 

The evolved three body state is 

U{t)\* GHZ ) = V 1+ (t)®V 2+ (t)(ac\++)+bd*\—)) 12 \+} 3 

+Vi_(t) ® V 2 _(t)(ad\++) - 6c* I— ))i 2 |+) 3 ( 2 - 31 ) 
= Ep±Vi ± (t)®V2±(t)\x±)u\±h, 
where p± are the same as in Eq. 23 and 

\ X+ ) = l/p+{ac\++) + bd*\ )), 

= l/p-(ad\++)-bc*\—)), 
are the entangled states for the joint system 1,2. The reduced state is then, 

Pxi(t) = EP±^(*) ® V 2± (t)\ X ±)(x±\Vi ± (ty ® V 2± (ty. (2.33) 

The evolved state is a statistical mixture of entangled states, which depending on the 
weights of the mixture may yield greater than zero entanglement for the joint system 1,2. 
Therefore £{pi 2 (t)) >£(pi 2 (0)). 

b. Triple states 
Let 

I**) = a|001) + 6|010) + fc|100) 

= RiR 2 (a\++) 12 R 3 \-} 3 + b\+-) 12 R 3 \+} 3 + k\++} 12 R 3 \+} 3 . 
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The initial tangle between qubits 1 and 2 is r(pi 2 (0)) = 4|6A;| 2 . After evolution, the state is, 
U(t)\V tr ) = V 1+ (t)®V 2+ (t)(ad*\++} + bc\+-} + kc\-+}) 12 \+}3 



+Vi_(t) <g> V 2 _{t){-ac*\++) + bd\+-) + M-+))i2|+)s- 



(2.35) 



where we have used the rotation R 3 defined in |277|. The reduced state is 

p 12 (t) = Y / mlV 1± (t) (g) K 2± (t)|0 ± )(0 ± |V 1± (t)t ® \/ 2± (t)t, (2.36) 

where 



m + = v l a r + l c r — 2|a| 2 |c| 2 , 
m_ = ^|a| 2 + |d| 2 -2|a| 2 |(i| 2 , 
and the normalized states for qubits 1 and 2 are 

= l/m+(ad*\++) + c&|+-) + cfc|-+)), 
= i/ m _(_ ac *|++)+tZ6|+-)+dA;|-+)). 



(2.37) 



(2.38) 



(2.39) 



In order to relate the entanglement between qubits 1 and 2 before and after evolution, 
we use the fact that the tangle is a convex function on the set of density matrices |5[] such 
that 

r(pi2(t)) < E™ 2 ± r(|0 ± » 

= m 2 + (A\cb\ 2 \ck\ 2 /m 2 + ) + m 2 _ (4\db\ 2 \dk\ 2 /m 2 + ) 
= r(pi 2 (0))(|c| 4 + (l-|c| 2 ) 2 ) 
< r(p 12 (0)). 
Therefore £{pn{t)) <£{pu(0)). 

It is natural to ask what can be said of the residual three tangle under this evolution. 
Coffman et.al. |J have constructed a set of three quantities (di, d 2 , d 3 ), invariant under 
permutation of three qubits, that can be used to quantify the residual tangle as ri 23 = 
4|g?i — 2d 2 + 4(i 3 |. In terms of a state expanded in the logical basis |£) = /ttnk\ijk), they 

ijk 

are: 

di — a ooo a m + a ooi a iio + a oio a ioi + a ioo a oii' 

d 2 = aooo a iiifloiiflioo + a ooofliiiflioiOoio + a oooOiiiaiioaooi .„ 

(2.40) 

+aonaioo«ioiaoio + ooiifliooflnoOooi + aioiaoioOnoOooi; 

^3 = OooO a 110 a 101 a 011 + a lll a 001 a 010 a 100- 

Mapping |+) — > |0), | — ) — > |1), it is clear from Eq. |2.4| that the residual tangle may 
either increase or decrease over time and will return to its initial value at times \a\t = kn/2 
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— * 

if = k/l for some integer pair k,l. The only invariant quantity above is o^. If 

one identifies the matrix elements as vertices of a Boolean cube then d 3 is a sum of 
two configurations each lying on the vertices of a tetrahedron ||. These configurations 
are composed of four states with all even or odd Hamming weight. For states with even 
or odd Hamming weight, d\ and d 2 are zero and from the definition of r 12 3 the residual 
tangle is conserved under the evolution from commuting Hamiltonians with value r 12 3 = 
16|aooo a iio a ioi a oii| or r i23 — 16|aiiiaooi a oio a ioo|- 



III. NONCOMMUTING INTERACTIONS 



In contrast to the case for commuting Hamiltonians, noncommuting pairwise interactions 
can increase or decrease the entanglement between the noninteracting pair given any class 
of initial states. To illustrate this we consider the Heisenberg interaction on the ID chain of 
qubits oriented (1,3,2) with equidistant spacing: 

<i,j> 

= ga x • 0-3 + ga 2 ■ a 3 ( 3 - 1 ) 
= H 13 + H 23 . 

There are three distinct eigenvalues of H namely (Eq — 0, E\ — —4:g, E2 = 2g) corre- 
sponding to eigenvectors antisymmetric under 1^2, symmetric under 1 <-> 2, and symmet- 
ric under all permutations respectively. Any state that is permutation symmetric, such as 
a tensor product of three identical states, is stationary meaning the entanglement between 
qubits 1 and 2 is invariant. However, the state = 1 00+) , for instance, has overlap both 
with eigenstates invariant under permutations of 1, 2 and 1, 2, 3. The reduced state of qubits 
1 and 2 can be easily be calculated by diagonalizing the Hamiltonian and finding the Krauss 
operators A = 3 (0\U(t)\+) 3 and A 1 = 3 {l\U(t)\+) 3 as, 

Pn(t) = A (t)|00)(00|A (t) t + A 1 (t)|00)(00|A 1 (t)t 

= (2/9(1 + cos(6^)) + 5/9)|00)(00| + 2/9(1 -cos(6^))|V' + )(V' + l ( 3 - 2 ) 
+V2/3ism(3gt)e- i39t \00)(ip + \ - V2/3i sin(3#t)e i3ff *|V> + ) (00|. 

The tangle of this state is r^pnif)) = 4/9 sin 3 (3gt). 

For the other classes of states we do not present the calculations but merely state results. 
For initial bipartite entanglement between 1 and 2, the tangle r(pi2(t)) can be either greater 
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or less than r(pi 2 (0)), while when the entanglement is initially shared between 1 and 3 
(or 2 and 3), r(pi 2 (t)) > i~(pi 2 (0)) = 0. The tangle is zero for generalized GHZ states 
because these states are permutation symmetric. The triple states can have r(pi 2 (t)) either 
be greater than or less than t(/? 12 (0)). It is notable that the residual tangle is invariant for 
the triple states (r 12 3 = 0) and the generalized GHZ states (r 123 = 1). This is a property 
of the Heisenberg interaction having permutation symmetric eigenstates corresponding to 
these two classes. Such symmetry persists even when the coupling constants differ but 
is broken for any anisotropy in the interaction. For more general states the residual tangle 
can either increase or decrease under evolution. 



IV. CONCLUSIONS 

We have investigated the preparation of entanglement between two qubits, 1 and 2, 
that interact pairwise with a third qubit 3 but not with each other and in principle could 
be quite far apart from one another. By explicit calculation we have found the change 
of entanglement between 1 and 2 for various classes of initial states when the interaction 
Hamiltonians commute. This is the case when the time order of the interactions is irrelevant. 
Under such evolution, the entanglement of formation between 1 and 2 is nonincreasing for 
all initial pure states except the generalized GHZ states where it can increase from zero. In 
contrast, noncommuting Hamiltonians can generally increase or decease the entanglement 
between the noninteracting pair of qubits. 

One may be able to prepare entangled states from initial product states by performing a 
measurement on the third qubit. This, in fact, is at the heart of the entangled chains proposal 
for quantum computing by Briegel and Raussendorf [|l!| wherein a neighborhood of qubits 
perhaps confined to a lattice are allowed to interact by mutually commuting o % z o\ interactions 
and arbitrary subsets can be prepared in entangled states by measurements on the remainder 
of the neighborhood. In that proposal, because of the nature of the commuting interaction, 
measurement and local addressing of the qubits is necessary to entangle distant pairs. Any 
entangling interaction can be mapping to any other two body entangling interaction by 
local unitary operations [13], however, if within a given system single qubit addressability is 



difficult this may not be possible. A focus of future research will be to study what kind of 
entangled states could be synthesized given finite measurement resolution over sets of qubits 
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interacting pairwise with each other. 
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